QUANTITATIVE EQUIDISTRIBUTION FOR THE SOLUTIONS OF 
A SYSTEM OF SPARSE POLYNOMIAL EQUATIONS 

CARLOS D'ANDREA, ANDRE GALLIGO, AND MARTIN SOMBRA 

Abstract. For a system of Laurent polynomials /i, . . . , / n £ C[xj 1 , . . . , x„] whose 
coefficients are not too big with respect to its facet resultants, we show that the 
solutions in the algebraic torus (C x ) n of the system of equations /i — ■■■ — f n — 0, 
are approximately equidistributed near the unit polycircle. This generalizes to the 
multivariate case, up to a exponent, a classical result due to Erdos and Turan on 
the distribution of the arguments of the roots of a univariate polynomial. 

We apply this result to bound the number of real roots of a system of Laurent 
i3 \ polynomials, and to study the asymptotic distribution of the roots of systems of 

Laurent polynomials over Z and of random systems of Laurent polynomials over C. 
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1. Introduction and statement of results 

A celebrated result due to Erdos and Turan says that, for a univariate polynomial 

over C whose middle coefficients are not too big with respect to the extremal ones, 

>-^ ■ the arguments of its roots are approximately equidistributed |ET50| . Combined with 

,^Zh I a recent result of Hughes and Nikeghbali |HN08| , this shows that the roots of such a 

C3 ■ polynomial cluster near the unit circle. 

We introduce some notation to precise this result. Let Z be an effective cycle of C x 
of dimension 0, that is, a formal finite sum 

Z = /Z m ^ 

with £ G C x and m^ G N. The degree of Z, denoted deg(Z), is defined as the sum of 
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its multiplicities m^. We assume that Z ^ or, equivalently, that deg(Z) > 1. 
For each — ir < a < [3 < ir, consider the cycle 

Za,f3 = Yl m ^' 

a<arg(£)</3 

where arg(£) denotes the argument of £. The angle discrepancy of Z is defined as 

deg(Z aj/3 ) j3 - a 



^ A ang (Z) 



sup 



deg(Z) 2vr 



-?r<o</3<7r 

For < e < 1, consider also the cycle 

z £ = Y2 m ^ 



l-e<|$|<(l-e)- 
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The radius discrepancy of Z with respect to e is defined as 

deg(Z) 

For a polynomial / G C[x] \ {0}, we denote by Z(f) the 0-dimensional effective cycle 
of C x defined by its roots and their corresponding multiplicities. We also set ||/|| S up = 
SUP| Z |=1 |/0)l- 

Theorem 1.1. Let f = ao + - • ■ + a ( ix d G C[x] with d > 1 and ao°d 7^ 0, and < e < 1. 

T/ien 



1 



A ang (Z(/))<c^ -log 



sup 



A rad (Z(/),e)< — log. 



sup 



wi/l c = y / 27r/K 
Zan's constant. 



2.5619 . . . , where k = ^ 



oo (-ir 



m=0 (2m+l) 2 



0.915965594. 



is CWa- 



The more interesting (and hardest) part is the bound for the angle discrepancy. 
The original Erdos-Turan result states [ET50| 



A ; , 



5 (Z(/))<16W-log 
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A few years after that paper, Ganelius |Gan54j replaced the £ 1 -norm ]T\- \a,j\ by the 

smaller quantity ||/|| S up an d improved the value of the constant to c < \j2i\ j k. On 
the other hand, Amoroso and Mignotte |AM96j showed that the optimal value of c 
cannot be smaller than y2- The bound for the radius discrepancy is due to Hughes 
and Nikeghbali |HN08| . 

In this text, we study the distribution of the solutions in the algebraic torus (C x ) n 
of a system of multivariate polynomial equations. For instance, consider the following 
system of bivariate polynomials: 



(1.2) f 1 = x\ i + Xl xY + x^ + l, h = x\ 2 x 2 - x\ 6 - x x x 2 + 1 eC[n,i 2 ], 

These are polynomials with a moderate degree and small integer coefficients. By direct 
computation, we can verify that the solutions in (C x ) 2 of the system of equations f\ = 
/2 = are aproximately equidistributed near the unit polycircle S 1 x 5 1 (Figure [1]). 
This example and others of the same kind suggest that Theorem 11.11 has an extension 
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Figure 1. Angle and radius distribution of the zeros of the system ([1.2 
to higher dimensions. 
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The study of the distribution of the solutions of a system of multivariate polynomial 
equations has been addressed from different perspectives. In the frame of his theory 
of fewnomials, Khovanskii has shown that one can estimate the distribution of the 
arguments of these solutions in terms of the number of monomials and the Newton 
polytopes of the input system, see for instance |Kho9H §3.13, Theorem 2]. There are 
also several interesting results by Shiffman, Zelditch and Bloom on the asymptotic 
distribution of the solutions of a random system of polynomial equations, see for 
instance [SZ041 IBS07J and the references therein. 

Our purpose in this text is to obtain an extension of Theorem 11.11 to systems of 
Laurent polynomials with a given support. For i = 1, . . . ,n, let Ai be a non-empty 
finite subset of Z n and Q { = conv(A) C R n its convex hull. Set D = MV(Q X , ...,Q n ) 
for the mixed volume of these lattice polytopes, and assume that D > 1. For each i, 
let fi be a Laurent polynomial with support contained in Ai, that is, 

fi= X^^a^eC^ 1 ,...,^ 1 ] 

with a.i :a G C and x a = x^ 1 ...x^ for each a = (ai,...,a n ) G Ai- Write / = 
(/i, . . . , f n ) for short. For each primitive vector v G Z n , consider the corresponding 
"facet resultant", denoted by Res.4« ...,A v (fii ■ ■ ■ > fn) £ C, of / in the direction of v 
(Definition 13. 4p . These resultants can differ from 1 only for a finite number of v's. We 
assume that all of them are nonzero. This condition holds for a generic choice of / in 
the space of coefficients and, by Bernstein-Kushnirenko-Khovanskii's theorem jBer75|. 
Theorem B] , it implies that all the solutions of /i = • • • = f n = are isolated and 
that their number, counted with multiplicities, is equal to D. 

For a vector w G S n ~ l in the unit sphere of W n , we denote by w^~ C M. n its 
orthogonal subspace and by tt w : W l — > w^- the orthogonal projection onto it. We 
denote by MV^± the mixed volume of convex bodies of w 1 - induced by the Euclidean 
measure on w and, for i = 1, . . . , n, we set 

D W: i = MV w i (vr^(Qi),. . . ,7C w (Qi-i),n w (Qi + i), . . . ,n w (Q n )). 

We then define the Erdos- Turdn size of / as 

(1.3) r,(f) = ±- sup log' :.=illMlsup \ 



uuj, ^ 6 \(V, W )\ 1 

D weS n~i VnjRes^,...,^(/f,...,/-)i^y 

where the second product is taken over all primitive vectors o£Z™ and (•, •) denotes 
the standard inner product of W 1 . This is a generalization to the multivariate case of 
the quantity \ log I "f" sup ) that appears in Theorem ll.il (Proposition 13.121) . 

a \V|aoad|/ 

We set Z(f) for the O-dimensional effective cycle of (C x ) n defined by the roots 
of / and their multiplicities. The angle and radius discrepancies of cycles of (C x ) n 
are the obvious generalization of those for the univariate case (see Definition 12.11) . 

The following is our main result in this text. 

Theorem 1.4. Let A\, . . . ,A n be family of non-empty finite subsets of I/ 1 such that 
MV(Qi, ...,Q n ) > 1 with Qi = conv(A)- Let fi,...,f n G Cfccf \ . . . , x± x ] with 
supp(/j) C Ai and such that Res.4« ...,A%{fi> ■ ■ ■ > fn) ¥" f or a ^ v £ Z n \ {0}. Then 

A ang (Z(/)) < 66n2"(18 + log + ( ?? (/)- 1 ))l^ 1 )r / (/)i 
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with log + (x) = log(max(l,x)) for x > 0. Also, for < e < 1, 

In 
A rad (Z(/), £ )<— 77(/). 

For n = 1, our bound for the angle discrepancy is worse than Erdos-Turan's, due 
to our approach. By analogy with the one-dimensional case, it is natural to ask if, in 
the setting of our result, a stronger inequality of the form 

holds, with c(n) > not depending on /. It would be interesting to settle this 
question. 

Theorem II .41 has several consequences. For instance, we can derive from it a bound 
for the number of positive real solutions of a system of polynomial equations, in terms 
of its Erdos-Turan size. For a cycle Z = J2t 7n 4^ OI " (C x ) n , set 

Z + = Y^ m t€- 

£e(R> )" 
Corollary 1.5. Let notations be as in Theorem \1.4\ Then 

deg(Z(/)+) <66n2"(18 + log + (7 ? (/)- 1 ))K n - 1 )7 ? (/)^ deg(Z(/)). 

We can also apply our result to study the asymptotic distribution of the roots 
of a sequence of systems of polynomials over Z with growing supports and whose 
coefficients are not too big. To be more precise, let Qi, i = l,...,n, be lattice 
polytopes in ~R n such that MV(Qi, . . . ,Q n ) > 1- For each d > 1 and i = 1, ... , n, 
consider the finite subset of Z n given by 

(1.6) Ai, d = dQinZ n . 

For a Laurent polynomial / € C[x 1 , ...jX^ 1 ], we denote by supp(/) its support, 
defined as the subset of Z n of the exponents in its monomial terms. We also set 

ll/llsup = SUp K | =lv .. iK | =1 \f(wt, . . . ,W n )\. 

For a nonzero cycle Z = ^£ m £^ of (C x ) n , we consider the discrete probability 
measure on (C x ) n defined by 

where o^ is the Dirac measure supported on the point £,. Let v be the measure 
on (C x ) n supported on (S l ) n and whose restriction to this polycircle coincides with 
its Haar measure of total mass 1. 

Recall that a sequence of measures {vd)d>i on (C x ) n converges weakly to v if, for 
every continuous function with compact support h: (C x ) n — > M., it holds 

lim / hdvd = / hdu. 

d->oo 7(C x ) n J(c x ) n 

If this is the case, we write lim<2->oc Vd = v - 

Theorem 1.7. For d > 1, let f d = (fi,d, ■ ■ ■ ■> fn,d) be a family of Laurent poly- 
nomials in 7L\x 1 ,...,x^] such that supp(/j d) C dQi, log ||/id|| S up = o(d), and 
Res ^>-.< d (/l> ■ • ■ . /n,d) ^ /or a// v G Z» \ {0}. Tften 

lim o" z(f ) = i/. 

d^oo (Jd> 
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This result admits a quantitive version, which we state in Proposition 14.21 Theo- 
rem 11.71 is related to Bilu's equidistribution theorem for the Galois orbit of algebraic 
points in (C x ) n of small height |Bil97j which, at least for n = 1, also admits quantita- 
tive versions |Pet05l IFRL06J . Indeed, when the cycles Z(f d ) are irreducible over Q or 
when, at least, they have a very large irreducible component, it is possible to deduce 
Theorem 11.71 from Bilu's theorem and the arithmetic Bezout theorem. 

We can also apply Theorem 11.41 to study the distribution of the roots of a ran- 
dom system of Laurent polynomials over C We will show that, under some general 
conditions and without assuming any independence or equidistribution condition on 
the coefficients of the system, these roots tend to cluster uniformly near the unit 
polycircle. 

Let us keep notations as above and set A. d = (^4i,d, ■ ■ ■ , -4n,<i) with Ai t d as in p. 61) . 
Each point of the projective space P(C d ) can be identified with a system f d = 
(fi,d, ■ ■ ■ , fn,d) of Laurent polynomials such that supp(/j^) C dQi, i = 1, . . . , n, mod- 
ulo a multiplicative scalar. The associated cycle Z(f d ) is well-defined, since it does 
not depend on this multiplicative scalar. 

Let \x d be the normalized Fubini-Study measure on P(C^ d ) of total mass 1, and g d a 
probability density function on P(C"^ d ) (see 3Hfor details). Let f d be a random system 
of Laurent polynomials with supp(fi >d ) C dQi, i = 1, . .. , n, distributed according to 
the probability law given by g d with respect to /i d . The expected zero density measure 
of f d is the measure on (C x ) n defined, for a Borel subset U, as 

E(Z(f d );g d )(U)= f deg(Z (/ d )\ v ) g d (f d ) dfi d , 

where Z{f d )\ v denotes the cycle T,^eV(f d ) nU m ^- 

Theorem 1.8. For d > 1, let g d be a probability density function on P(C d ) with 
respect to the measure fj, d , and f d = (/i,d, . . . , f n ,d) a- random system of Laurent 
polynomials with supp(/j i d) C dQi, i = 1, . . . , n, distributed according to the probability 
law given by g d . Assume that the sequence (gd)d>i is uniformly bounded. Then 

E(Z(f d );g d ) = 
d^ood n M\{Q 1 ,...,Q n ) 

As an application, consider a random system of Laurent polynomials f d with 
supp(/i d ) C dQi whose coefficients are independent complex Gaussian random vari- 
ables with mean and variance 1. The random cycle Z(f d ) might be described by 
the uniform distribution on P(C^) (see Examp le 1 4 . 1 4 1 f or details). Then, Theorem II .81 
implies that the roots of f d converge weakly to the equidistribution on (S 1 ) 71 , and we 
recover in this way a result of Bloom and Shiffman [BS07, Example 3.5]. 

Our strategy for proving Theorem 11.41 consists in reducing to the univariate case. 
In 32 we consider the problem of studying the angle and radius discrepancies of an 
arbitrary O-dimensional effective cycle Z in (C x ) n in terms of the angle and radius dis- 
crepancies of its direct images under all monomial projections of (C x ) n onto C x . By 
applying a tomography process based on Fourier analysis, we show that the distribu- 
tion of Z can be controlled in terms of the distribution of its projections (Theorem [23]) . 
In 33 we consider cycles defined by a system of Laurent polynomials with given sup- 
port and we compute their direct image under monomial projections, in terms of 
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sparse resultants. Theorem II .41 then follows by applying the univariate Erdos-Turan's 
theorem combined with Theorem \2.'d\ and the basic properties of the sparse resultant. 
In § HI we study the asymptotic distribution of the roots of a sequence of systems of 
Laurent polynomials over Z and of random systems of Laurent polynomials over C. In 
both situations, the key step consists in bounding from below the size of the relevant 
facet resultants. In the case of systems over Z, this is trivial since these facet resultants 
are nonzero integer numbers. In the case of random systems over C, the result follows 
from an estimate of the volume of a tube around an algebraic variety due to Beltran 
andPardo [BP07J . 

Acknowledgments. We thank Carlos Beltran and Michael Shub for useful discus- 
sions and pointers to the literature. Part of this work was done while the authors met 
at Universitat de Barcelona and Universite de Nice-Sophia Antipolis. We thank these 
institutions for their hospitality 

2. Angle and radius distribution in the multivariate case 

In this section, we show that the angle and radius discrepancies of an effective 0- 
dimensional cycle in the torus (<C X ) n can be bounded in terms of the angle discrepancy 
of its image under monomial maps from (C x ) n to C x . 

Let Z be a nonzero effective 0-dimensional cycle of (C x ) n , which we write as a 
finite sum 

with 7?i£ G N and £ G (C x ) n . The support of Z is the finite subset of (C x ) n defined 
as \Z\ = {£ | m£ > 1}, and the degree of Z is the positive number deg(Z) = £\ m^. 

Definition 2.1. Let Z = £\ m^$ be a nonzero effective 0-dimensional cycle of (C x ) n . 
For each a = (a±, . . . , a n ) and f3 = (/5i, . . . , fi n ) with — i: < ctj < /3j < it, j = 1, . . . , n, 
consider the cycle 

-aj<axg(£j)<Pj 

The angle discrepancy of Z is defined as 



Aang(-Z) = sup 

a,/3 



deg(Z Q , iy a) -fj /3j - a 



n 



deg(Z) 1_1 2t 



Let < e < 1 and consider also the cycle 

i- £ <|^|<(i- E )-i 

where £j is the j-th coordinate of $,. The radius discrepancy of Z with respect to e is 
defined as 

deg(Z £ ) 



A rad (Z,e) = l 



deg(Z) • 



We have < A ang (Z) < 1. Observe also that < A rad (Z, e) < 1 and A rad (Z, e) = 
for all e if and only if \Z\ C (S 1 ) 71 . 
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For a lattice point a = (oi,. . . ,a n ) G 7L n we denote by x°"- (C x ) n — > C x the 
associated character, defined as x a (£,) = £1 1 ■ ■ ■ 6n n f° r £ e (C x ) n . The direct image 
of Z under \ a is the cycle of C x given by 

i 

We also set 



(2.2) 



9(Z) 



AaDg( *\ (Z)) , p(Z >£ ) = f; A rad ( X ^(Z), £ ), 



aeZ n \{0} 



a 



i=i 



where e^ denotes the j-th vector in the standard basis of Z n , and ||<x|| 2 is the Euclidean 
norm of the vector a G Z n . We have < 6{Z) < 1 and < p(Z, e) < 1. 

Theorem 2.3. Let Z be a nonzero effective 0- dimensional cycle of (C x ) n . Then 

A ang (Z) < 22n(|V(9 - log^)))^" 1 ' 6{Z)l 

and, for < e < 1 , 

A rad (Z,e)<p(Z, £ ). 

The rest of this section is devoted to the proof of this result. For u,v £ W 1 we 
write {u,v) = X^?=i u j v j f° r their standard inner product, and for £ G (C x ) n we set 
arg(£) = (arg(£i), . . . ,arg(£ n )) G (-vr,7r] n . 

Lemma 2.4. Let Z be a nonzero effective 0- dimensional cycle of (C x ) n and a G 
Z n \{0}. T/ien 



deg(Z)^ * 



< 27rA ang ( X »(Z)). 



Note that for a = we get ^^E— E^ ™,£ e'* '^®* = 1. 

Proof. Set L>o = #1-^1 and L> = deg(Z) for short. Let — it < v^ < ir, k = 1, . . . , £>o> 
denote the inner products (a, arg(£)) modulo 2tt for the different points £ in the 
support of Z, and let m& denote their corresponding multiplicity. We suppose that 
these numbers are arranged in increasing order, that is, v\ < ■ ■ ■ < vd - 
For — 7r < v < it set 

JV(l/) = ^2 m k- 

We have 



/7T A) „j, fc+1 / ^^ 

-7T fc=1 Ju k \ l<k 



mi I e lu dv 



Vk + l 



« + E 



m fc e 



Wk 



E(E 

where we have set ^d +i = 7r - On the other hand, an easy calculation shows that 



(2.6) 



e dv = \. 

.„ 2vr 



8 D'ANDREA, GALLIGO, AND SOMBRA 

Combining (|2.5|) and (|2.6|) . we deduce that 

D a 

D 



m ^ e i(a,arg(0> = j-^m k e 1Uk 
fc=l 



, A) 



,V-Mr_JVp. ^ 



2n 



D 



Hence 



^E-^ 1<a 



arg(^)> 



< 



U + TT N{v) 



du 



2vr D 

is + tt deg(x?(Z) ,„) 



2vr deg(x?(Z)) 



aV < 2vrA ang (x^(Z)), 



which concludes the proof 



□ 



Let a, (3, r G R such that a < j3 and r > 0. We consider the function /i 
defined, for x G R, by 



K,/3,t(x) 



if x < a — r, 
ff( x ~" +T ) if a — r < x < a, 

1 if a < x < (3, 
g(£±rzl) if/3<x</3 + T, 
if /3 + r < x, 



with g(x) = — 2x 3 + 3x 2 . Lemma 12.71 below shows that h a r is an approximation of 
the characteristic function of the interval [a,/3]. 

For m G N, we denote by C m (R) the space of functions /: R — > R having m con- 
tinuous derivatives. 

Lemma 2.7. Let a,/3,T € R swc/i i/iai a < (3 and r > 0. TTiera 

(1) h aAr G C^R); 

(2) h at p, T (x) = 1 for x £ [a, 13], h a ,p, T {x) = for x £ (-oo,a - r] U [/? + r,oo), 
and < h a fi )T (x) < 1 /or aZZ x G R; 

(3) X!^, ha 7 )T dx = (3 -a + r and, moreover, /"_ T h a> p >T dx = J^ +T h a ^ r dx = § ; 

(4) /^ l^, T |cfa = 2; 

(5) /i'^ a r is absolutely integrable and ff^ \h'^g T \dx = ~. 

Proof. By a direct calculation, we verify that the function g satisfies the following 
properties: 

• g(x) > Ofor allxG [0,1]; 

• 5(0) = </(0) = 0, g(l) = 1, e/(l) = 0; 

• Io9dx = i; 

• Jo l^'l d2; = !; 

• Jo \a"\dx = 3. 

The claim follows easily from these properties and the definition of h a g T . D 

Suppose furthermore that f3 — a + 2r < 2ir. The support of h aj p jT is then contained 
in an interval of length bounded by 2ir, and so this function can be regarded as a 
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function on R/27rZ. Consider its Fourier series 



h a ,p, T (x) ~ /^c a e h 



a6Z 



with c a = i f^ n h aj jT (x) e mx cfa. 



Lemma 2.8. Let a < )3 and r > such that /3 — a + 2t < 2-/r. T7te7J X^aez Ca emx 
converges absolutely and uniformly on R/27rZ to h a ^ :T . Moreover, cq = 2 " T , and, 
/or a ^ 0, 



(2.9) 



; 1 3-1 

-7ra ttt a z ) 



Proof. Lemma I2.7I| T]) implies that the series 2~^aez c « eia:C converges absolutely and 
uniformly on K/27rZ to the function /i Qj/ g ]T . The computation of cq follows from 
Lemma l2T7l![ 3|) . Integrating by parts, we deduce for a £ Z \ {0} that 



-2-7ria 



h'(x)e- iax dx 



27r(-ia) 2 ./_„ 



< j/3iT (x)e- i -^. 



Hence, IcJ < rk- f 71 ' l/i' « Jdx and also IcJ < t^-t f 71 ' \h" aJdx. Then (12.91) follows 
by bounding these integrals with Lemma [22U3HSD • D 

Next, we apply Fourier analysis to control the angle discrepancy of Z in terms of 
the angle discrepancy of its direct image under monomial projections. 

Lemma 2.10. Let n > 2, q £ Z>i, a = («i,...,a n ) and (3 = (/3i, . . . ,/3 n ) ureia 
ctj, f3j G M suca i/iai — 7r < ay < /3j < 7r and /3j — a,- H — < 27r. Taen 



deg(Z Qj/3 ) 



deg(Z) 



n 



& - a. 



2tt 



< 2?i6»(Z) 



3ro 

2vrg 



in+3 



V3 



+ n^^-o2(9 + log(g)) n -^(Z). 



7T 



Proo/. Set r = | and h a ,p, T {v) = I\]=ih aj ,^A u j) for v = ( u l,- ■ ■ , u n) G 
also D = deg(Z) and 

deg(Z aii3 ) 1 v^ , / ,>» 

— ^ ^ /^"^./^(argio) 



Set 



i=i 



2tt 



— ^ m^/ laAr (arg(^)) - JJ 

n 



ft ~~ a j + T TT ft ~~ a i 



2tt 



n 



2tt 



We will bound each of these quantities. For Si, we consider the subset of M. n given 
by /c,,/3,t = E["=i[«j - r > ft' + T l \ Ylj=i[^j,/3j]- Then 



1 E 

D *-" 



m 



C^a,/3,r(arg(^)) 



arg(^)e/ a , /3iT 



10 



D 'ANDREA, GALLIGO, AND SOMBRA 



For each £ such that arg(£) G IocBti there is 1 < j < n such that either otj — r < 
arg(^j) < Qj or j3j < arg(£j) < j3j + r. Since < h a ^^ T (v) < 1 for all v, we have 
that Si is bounded above by 

1 n 

n Yl ( de g(X* J ( Z )«i-r,a,) + deg(x* J (^) / 3 jI/ 3 J +r))- 



i=l 

By using the definition of A ang and Lemma 12.41 we get 

ideg( X ?(Z) a ,_ T , a .) < A ang ( X ?(Z)) + ^ < 0(Z) + ^ = 0(Z) + JL 
and a similar bound holds for deg(x* J (-^)/3 ,/3 + T )- Hence, 



(2.11) 



Si < 2n(9(Z) 



■/) 



7Tf/ 



Now we turn to S2. Due to the conditions imposed on r, we can regard h a ^ )T as 
a function on M Tl /27rZ™ ~ (— 7r,7r] n . Let Saez nCae be ^ s multivariate Fourier 

series. For j = 1, . . . , n, we denote with X^a-ez c J> a j eiajUj the Fourier series of h aj ^ jtT . 
Then, for each a = (01, . . . , a n ) £ Z n , 



n- 



-j,aj • 



In particular, cq = TT?— 1 J -T J T - The Fourier series of each h a . g r converges abso- 
lutely to this function, and so the same holds for the Fourier series of /i Q „3, r - Hence, 

for 1/ e (— 7r,7r] n . Applying Lemma 12.41 we obtain 



(2.12) S s 



^E^E <=.«' <wM,) sEw^E 



€ «7^0 






mt e 



i(a,arg(£)} 



<2TT8(Z)Y J \ca\\\a\\l <2 n+1 Tr9(Z)J2 (Y,^)\ c ° 

a^O a>0 ^ s=l 



For s = 1, 

(2.13) 



^\^T|Ca| = ( J^ V^I|C1, 01 |) n ( ^ |c i:a .| j 
a>0 ^ai>0 ' 7=2 ^a,>0 ' 
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Using the bounds in (|2.9p we get, for j = 2, . . . , n, 

3q 1 oo „ 



^ic w i<| Cj , | + ^ — + J2 



"1 a 3 n Yi, 7rr °? 

a,j=l J cij=Aq+l J 

fc-aj+T ( 1 f 3q dx\ 3 f°° dx 

2vr 1 , /3\ 3 r 

<^- + l + -log - +— - 



(2.14) <-(9 + logg). 

7T 

Similarly, we now bound 

3? -, oo _ 
V \/aT|ci, ai | < V — — + V — 

ai>0 ai = l v ai=3ij+l v 

^ 1 /" 3g _I , 3 f°° _3 , 

< — / x 2 dx -| / x idx 

* Jo KT J 3q 

4^/3 i 

(2.15) < -^—q*. 

TT 

It follows from (T2J2D . (I2JD and (I2J51) that 

on+3 /q 

(2.16) S 2 < „f_jV(9 + log^))"- 1 0(Z). 

Next we consider S3. Set 0(£) = JT? 1 3 2n ^ or * £ ^- There exists < to < 1 
such that S 3 = |0(1) - 0(O)| = \(f>'(t )\ and so 

S 3 < SUp 0(*o) <2^9~11 5 - nT To ^n = 0~- 

Finally, we collect (|2.1ip . (J2.16P and the above inequality to get 
deg(Z at/3 ) -A- fa - aj 



n 



deg(Z) 1 = 1 2vr 



<Si + S 2 + S a 



<2n0(Z) + - + n^^M (9 + log(g))"- 1 9{Z) + -£- 
=2n6(Z) + \-+ rS^qHv + log^))"" 1 9{Z). 

ITXq TT n i 

n 

Proof of Theorem \2.3[ For the radius discrepancy, we have 



\Ze\ 



n .. 

n{€€iziii- e <fei< I - i }. 
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By taking complements in this equality and considering the corresponding multiplic- 
ities, we deduce that 



v. 



deg(Z) - deg(Z £ ) < ^deg(x7(Z)) - deg( X 7(Z) £ ). 
i=i 

Hence, A rad (Z,e) < Y!j=i ^radixt 3 (Z , e)) = p(Z,e), as stated. 

We now consider the bound for the angle discrepancy. For n = l, A ang (Z) < 6{Z) 
and so the bound in the claim is trivial. Hence, we suppose that n > 2. Put then 
C(Z) = (9 - log(0(Z)))f ( n_1 ) 6(Z) 2 / 3 G M> for short and set 



9 §(n-l) 



Suppose also that q > 1. Then 



g^™- 1 ) 9a( n_1 ) 1 

<q< .._. < 



(2 - 17) 2C(Z) ^ C(^) "W 

Let a = (ai, . . . , a n ) and /3 = (/3i, . . . , /3 n ) with — it < aij < (3j < it. Consider first 
the case when fy — aj < n. In particular, f3j — ccj -\ — < 2ir. Applying Lemma 12.101 

and the inequalities (|2.17|) . we deduce that the quantity 
bounded above by 



dcg(Zq./3) 

deg(Z) 



in ft-ay 



In 



az) 



3ra 2C(Z) 2 n+3 



V3/9§^ 



'gl(n-i) ' 2 7rg|(n-i) ' " tt"" 1 V C( z ) 
Since n > 2, this quantity can be bounded by 



llj = l 2vr 

(9-log(^(Z)))^^(Z). 



is 



n 



9 §(n-l) 



+ 



)§("-!) 



+ 



7T 



7T 



ra-1 



<n 1 + 



C(3) 

2 3 v / 3tt/2- 



9^ 



7T 



C(Z) < 22 n 



«Z), 



A 2 3 x/3tt < 21 
93 



as it can be easily verified that — |— < | and 

If q = 0, then >,^ < 1, which implies that A ang (Z) < 1 < (,{Z). Hence, in the 
case when — n < oij < f3j < tt for all j, it holds 



(2.18) 



deg(Z 



a,0) 



deg(Z) 



IT 

3 = 1 



Pi ~ a 3 



2tt 



<22n(-J ((Z). 



Now, if (3j — aj > tt for a j, we subdivide each of those intervals (atj,f3j] into two 
subintervals of length < tt. We can then decompose Z a ^ as the sum of at most 2 n 
cycles of the form Z V0(T)l , la where the jth coordinate of vq^ (respectively Vi l(r ) is 

either aj (respectively " J 2 J ) or Qj 2 J (respectively /3j). Also, we can expand the 
product 

Pj-aij 



n 



2tt 
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as the sum of the volumes of sets n^liX *^ . 



l Va ' a 'i . 1/1 '' r 'J l From here, we easily get that 



deg(Z Q]/3 ) 



deg(Z) 



3=1 



2tt 



< 



£ 



deg(Z, 



^0,ff>"l,<i/ 



deg(Z) 



n 

i=i 



v\ 



<?.j 



U 0,a,j 



2tt 



and the claim follows applying the bound ()2.18p . which has to be multiplied by 2 n . 
Altogether, we get 



□ 



A ang (Z)<22nl-\ C(Z), 
which concludes the proof. 

3. Bounds for the discrepancy in terms of sparse resultants 

In this section, we consider cycles defined by a system of Laurent polynomials with 
given support. We first compute their direct image under monomial projections, in 
terms of sparse resultants. Then, we will derive Theorem 11.41 from the univariate 
Erdos-Turan's theorem and the results in the previous section. 

We recall the definition of the sparse resultant, as presented in |DS12j . This def- 
inition differs from the classical one which can be found, for instance, in [GKZ94. 
ICLO05[IFS93] . 

Let M ~ Z n be a lattice of rank n and N = M v = Hom(M, Z) its dual lattice. 
Set M R = M ® R and 2% = N <S> R. The pairing between b G M R and v E N R will 
be denoted by (6, v). Let Ao, . . . ,A n be a family of n + 1 non-empty finite subsets of 
M and put .4 = (Ao, ..., A n ). Let Ui = {«i, }aeA , » = 0, . . . , n, be a group of #A 
variables and set u = {uq, . . . , u n }. For each i, let i^ be the general polynomial with 
support Ai, defined as 

(3.1) Fi= J2 Ui, a x a e C[u][M), 

where C[u][M] ~ C[tt][x 1 , . . . , x^ 1 ] is the group algebra. 

Let Tm = Hom(M, C x ) ~ (C x ) n be the algebraic torus associated to M, and set 



W A ={(x,u)eT M x f[ 

i=0 

This is an irreducible subvariety defined over 

n 

^:T M x]]>(C A ) 



At 



j=0 



Fi( Ui ,x) = 0y 
of codimension n + 1. Let 
[JP(C" 

i=0 



-A, 



be the projection onto the second factor. The direct image W_a_ under this map is the 
Weil divisor of nr=o IP ( C ' 4i ) S iven b y 



v*(W A ) 



deg(7r| w ^ A )7r(W^) if codim(7r(W^)) 




1. 



if codim(vr(IF^)) > 2. 

The sparse resultant associated to A., denoted Res^j., is defined as any primitive equa- 
tion in Z[u] of this Weil divisor. It is well-defined up to a sign. When tt*(W a ) = 0, 
we set Res„4 = 1 for convenience. 

According to the definition above, sparse resultants are not necessarily irreducible. 
If, as in |DS12j . we denote by Elim.4. what is classically called the sparse resultant (for 
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instance in [GKZ941 ICLU051 IFSSSj ). then Res.4 + 1 if and only if Elim.4 / 1 and, if 
this is the case, 

Res.4 = =t Elim^ A . 

We recall the basic properties of the sparse resultant that we will need in the 
sequel. For simplicity, we will restrict to the case when M = N = Z n , so that each 
Ai is a finite subset of 7L n . Consider the convex hull Qi = conv(.Aj), which is a lattice 
polytope of W 1 . The mixed volume of convex bodies P\, . . . ,P n of W 1 is defined as 

n 

MV(P 1 ,...,P n ) = J2(- 1 ) n ~ J E vol(P ll+ --- + P l3 ). 

j=l l<ii<---<ij<n 

Write C[x ] = C[x 1 , . . . , £„] for short. The height of a Laurent polynomial / = 
SaeZ" a a^ a G Cfar 111 ] is defined as 

h(/) = log(max|a a |). 

Proposition 3.2. Let notations be as above. Then, the sparse resultant is a multiho- 
mogeneous polynomial of partial degrees 

deg Ui (Res.4.) = MV(Q , • • • ,Qi-i,Qi+l, ■ ■ ■ ,Qn) 

and height bounded by 

n 
h(Res^) < Y, MV(Q , • • • , Qi-i, Qi+i, ■ ■ ■ , Qn) log(#A). 

i=0 

Proof. The formula for the partial degrees is classical, see for instance [GKZ94]. In the 
present form, it can be found in [DS12J . The bound for the height is given by |Som04l 
Theorem 1] for the case when the resultant depends on all the groups of variables 
ito, ••• , u n . The general case can be found in [DS12| . □ 

For a family of Laurent polynomials fi G Cfa; 1 * 11 ] with supp(/j) C Ai, i = 0, . . . , n, 
we write 

Res^.(/ ,...,/„) 

for the evaluation of the resultant at their coefficients. The following is the multiplica- 
tivity formula for sparse resultants. 

Proposition 3.3. Let A. = (Aq, ■ ■ ■ ,A n ) and A! = (A' ,...,A' n ) be two families 
of non-empty finite subsets of Z n and set A + A' = (Aq + A' , . . . , A n + A' n ). Let 
fi, f[ £ C[sc ], i = 0, . . . , n, be Laurent polynomials with support contained in Ai and 
in A\, respectively. Then 

Res A+A'(fofo, ■■■, fnf'n) = Resx(/o, ■ • • , fn) ■ Res^/o, • • • , f' n ). 

Proof. The case when the resultant depends on all the groups of variables uq, . . . , u n 
is done in [PS93( Proposition 7.1]. The general case can be found in |DS12j . □ 

The support function of a compact subset P C M n is the function hp : W l — > M. 
defined, for v € W 1 , as 

hp(v) = inf (a,v). 
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Let Bbea finite subset of Z n and / = Yl,b&B ftb% b a Laurent polynomial with support 
contained in B. For v G M. n , we set 

B v = {beB\ (b,v) = h conv(B) (v)}, f v =J2 ^x b . 

beB-" 
A vector v G Z n is primitive if the gcd of its coordinates is equal to 1. 

Definition 3.4. Let Ai, . . . ,A n C Z" be a family of n non-empty finite subsets, 
v G Z n \{0}, and v 1 - C M n the orthogonal subspace. Then ZTlu 1 is a lattice of rank 
n — 1 and, for % = 1, . . . , n, there exists b^ v G Z n such that Af — bi jV C Z n n v -1 ". The 
/acei resultant of «4i, . . . , «4. n in the direction of v, denoted Res,4« ,4«, is defined as 
the resultant of the family of finite subsets Af — b lv cPfl v 1 - . 

Let fi G C[x x , ...jX^ 1 ], i = l,...,n, with supp(/j) C A{. For each i, write 
/? = x bi < v gi tV for a Laurent polynomial ^ )t) G C[Z n fl f -1 ] ~ C[yf 1 , . . . , y^-i] with 
supp(5Tj )W ) C *4" — bi tV . The expression 

Res^,...,^(/r,...,/^) 

is defined as the evaluation of this facet resultant at the coefficients of the gi iV . These 
constructions are independent of the choice of the vectors bi v . 

We have that Res^,...,^ 7^ 1 only if v is an inner normal to a facet of the Minkowski 
sum X^?=i Qi- ^ n particular, the number of non-trivial facet resultants of the family 
A±, . . . ,A n is finite. 

With notations as in Definition 13.41 write / = (/1, . . . , f n ) for short. We denote by 
V(/)o C (C x ) the set of isolated solutions in the torus of the system of equations 
fi = ■ ■ ■ = f n = and we set 

Z(f)= Yl m ult(£|/)£ 

€6V(/)o 

for the associated O-dimensional cycle, where mult(£|/) denotes the intersection mul- 
tiplicity of / at a point £. For /o G Cfa; 1 ' 11 ], we set 

/ (z(/))=n/o(o muit( ^ i/) - 

The following result is known as the Poisson formula for sparse resultants. 

Proposition 3.5. Let A = {Aq, . . . ,A n ) be a family of non-empty finite subsets of 
I/ 1 and fi G C[x ] a Laurent polynomial with supp(/j) C Ai, i = 0, . . . ,n. Suppose 
that Res^,...,^(/f , ...,/») ^ /or all v€Z n \ {0}. Then 



Res A (f J 1 ,...J n )=±(]jRe SA v_ A v(f?,...J%)- h ^ v Af (Z(f 1 ,...,f r , 



))• 



the product being over all primitive elements v G Z n . 



Proof The case when the resultant depends on all the groups of variables uq, . . . , u n 
can be derived from JPS93] (see also |Min03| ). The general case can be found in [DS12] . 

D 

From now on, we fix a family of non-empty finite subsets Ai, ■ ■ ■ , A n of Z n such that 
MV(Qi, . . . , Q n ) > 1, with Qi = conv(^4j). We consider also a family of Laurent poly- 
nomials fi, . . . , f n G C[s ±1 ] with supp(/j) C Ai, such that Res^y,...,^ (/f , ...,/£)/ 
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for all v € "L n \ {0}. By Bernstein-Kushnirenko-Khovanskii's theorem |Ber75|, Theo- 
rem B], 

deg(Z(/))=MV(Qi,...,Q„)>l. 

Consider the projection 7r ei : M n — > W 1 ^ 1 given by 7r ei [x\,X2, ■ ■ ■ , x n ) = (x2, ■ ■ ■ ,x n ). 
If we regard each fa as a Laurent polynomial in the group of variables x' := (x2, ... ,x n ) 
with coefficients in the ring C[xi ±:L ], its support with respect to x' is contained in the 
finite subset Tv ei (Ai) of R n_1 . Set then 

R(f) = Res^ i(A)v ..^ i( ^ n) (f 1 {x 1 ,x'),...,f n (x 1 ,x')) £C[if] 

for the evaluation of the resultant at these coefficients in C[x 1 ]. 

Recall that the sup-norm of a Laurent polynomial / £ Cfa^ 1 ] is defined as ||/|| S up = 
sup^ 6 (5i)n |/(xo)|. In general, it holds that 

(3.6) h(/)<log||/|| sup 



This is a consequence of Cauchy's formula for the coefficients of the Taylor expansion of 
a holomorphic function on C^ (see for instance |Som041 page 1255]). For a hyperplane 
H C R n , we denote by MV# the mixed volume associated to the Euclidean volume 
on H, defined for families of n — 1 convex bodies of H. 

The following result gives a bound for the sup-norm of R(f). Its proof is a variant 
of that for |Som041 Lemma 1.3]. 

Lemma 3.7. Let notations be as above. Then 

n 

log[|i2(/)|U P < ^MV^dvr^^l^Ologll/.llsup. 

3=1 

Proof. Let k > 1. Proposition 13.31 implies that 

(3.8) R(f) kn = Res k ^ i{Al) ^ Mei{An) (f 1 (x l ,x')\...J n (x 1 ,x f ) k ), 

where kn ei (Aj) denotes the pointwise sum of k copies of n ei (Aj). For short, set 
R k = RsSivK ei (Ai),...MeM n ) and f k = (fi,...,fn), so that the identity above can 
be rewritten as R(f) k = Rk{f )■ By Proposition 13.21 the partial degrees of this 
resultant are given by 

deg u .(R k ) = MV e x({fc7r cl (Q*)}^i) = P^MV,^^^,)}^), 

where Uj is a group of j^kAj variables, for j = 1, ... ,n. In particular, the logarithm 
of its number of monomials is bounded above as 

log(# supp(i? fc )) < log I [ [ I J 

n 

< £ deg u . (ii: fe ) log(#£^ + 1) = Oik^ 1 log(fc + 1)), 
i=i 

since j^kAj < (k + l) Cin for a constant ci independent of fe. By Proposition 13.21 the 
height of this resultant is bounded above by 

n 

HRk) < J2 de S Uj (Rk) log(#kAj) = O^' 1 log(k + 1)). 

3=1 
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Let w\ G S l . Using (|3.8|) . (|3.6p . the previous bounds, and the fact that ||/ 3 fc || S up = 
1 1 /ill sup> we deduce that 

^log||i?(/)|| sup = log||^(/ fc )|| sup 

n 

< h(R k ) + ^deg Uj ( J R J )log||/ J fe || sup + log(#supp( J R fc )) 
i=i 

n 

= ^(^MV e x({7r ei (Q,)}^)log||/,|| sup )+0(A ; n - 1 log(fc + l)). 
i=i 

The result then follows by diving both members of this inequality by k n and letting 

k ->■ oo. D 

For a G Z™ \ {0} and z an additional variable, set 

E a {f) = Res {0)O}) ^ lr .. :Ai (2: - x a , /i, . . . , /„) G C[z]. 

Due to the Poisson formula for sparse resultants given in Proposition [33J we have that 
Z{E a (f)) = x*(Z(f)), and so E a (f) can be regarded as an elimination polynomial 
for the cycle Z(f) with respect to the monomial projection x° '■ 

By [PS12] . there exists m G Z such that £ ei (/)(xi) = x^Rtf) (see also [CLOPS! 
Proposition 5.15, Chapter 3] for a proof for the resultant of dense polynomials). Hence, 
Lemma [3 .71 can be regarded as a bound for the sup- norm of the elimination polynomial 
E ei (/). Our next step is to extend this result to an arbitrary a. Recall that 7r a : M. n — > 
a denotes the orthogonal projection onto the hyperplane a C R n . 

Lemma 3.9. Let notation be as above. Then 

n 

log||£a(/)||su P < ||a|| 2 ^MV a x({7r a (Q^)}^)log||/ i || sup . 

i=i 

Proof. Consider first the case when o 6 Z" is primitive. The quotient Z n /aZ is 
torsion- free and so a can be completed to a basis of Z n . Equivalently, there is an 
invertible matrix A G SL n (Z) with first row a. Set a, a,2, ■ ■ ■ , a n and b\,b2, ■ ■ ■ ,b n for 
the rows of A and of A -1 , respectively. There is a commutative diagram 




where ipA and ip^-i are the monomial isomorphisms given byxH- x = (x ai , . . . , x an ) 
and x I—?- x = (x bl , . . . , x bn ), respectively. Let y = (y±, . . . , y n ) denote the coordi- 
nates of the torus below. For £ = 1, . . . , n, set 

ff(y) = ip\-,h = h(y b \- ■ -,y bn ) e Cfor"], 

so that (ifA)*Z(f) = Z(f ). Hence, E a (f) = E*a e and so Lemma [3T71 implies that 



(3.10) log||£ a (/)|U P < ^MV e x({7r ei (N(//))}^)log||/- 

i=i 



An 

I Slip 5 
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where N(//) is the Newton^ polytope of fo. We have vr ei (N(//)) = B(N(f e )) C 
B(Q e ) for the linear map B: W 1 -> M™" 1 given by B{x) = ({x,b 2 ), . . . ,(x,b n )). 
Let {v 2 , ■ ■ ■ ,v n } be an orthonormal basis of a and consider a second commutative 
diagram 

C 

a- 1 - s- a- 1 

u 

where C is the linear map defined by C(y2, • • • , y n ) = 2/2^2 + • • • + y n v n - It is easy to 
verify that U E GL(o ) is uniquely determined by 7r a , B and C. Since C maps the 
canonical basis of W 1 ^ 1 into an orthonormal basis of a , it is an isometry between 
these two spaces. On the other hand, a straightforward computation shows that 

n 

u ( v i) = ^2( v j, b k)v k = bj, j = 2, . . . ,n. 

We note that b 2 , . . . , b n is a basis of the Z-module a^~ n Z n . The Brill-Gordan for- 
mula |Bou70l Chapitre 3, § 11, Proposition 15]) implies that vol(a J -/(a- L nZ n )) = 1 1 ct 1 1 2 - 
Hence, 

I det(C/)| = vol(a ± /(a ± n Z n )) = \\a\\ 2 . 
Standard properties of the mixed volume imply that, for j = 1, . . . , n, 

MV e x({fl(Q/)}^ i )=MV x({Cofl(Q / )}^ i ) 

= MV a x({c7o 7 r a (Q^)}^ i ) = \\a\\ 2 MV a ±({ir a (Q t )Wj)- 

In addition, ifA gives an automorphism of (S 1 )™ and so ||//|| S up = ||/H|sup- We 
conclude that, when a is primitive, 

n 
(3.11) log||£ a (/)||sup < l|a||2^MV a x({7r a (Q^}^ i )log||/ i || sup . 

i=i 

Now let a £ Z n \ {0} be any vector. Choose a primitive a 1 G Z n and m € Z>i such 
that a = ma'. Using Proposition I3.3[ we deduce that 

E a (f)(z) = ± J] E a ,(f)(uz) 

where /U m denotes the set of m-th roots of 1. Hence 1 1 ct 1 1 2 = tt*. 1 1 c*-' 1 1 2 , vr a = 7r a / and 

log ||-E7 (/)||sup < mlog ||-E'a'(/)llsup- The result follows from the bound (|3.1ip applied 
to a'. □ 

Proof of Theorem \1.4\ Let a £ Z n \ {0}. Applying Proposition 13.51 with fo = z — x a 
and fo = z — x~ a , we get that the product of the leading and the constant coefficients 
of E a (f) is equal to 



±l[Res A1 ,... Al ,(f ?,..., a^\. 
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Recall also that Z(E a (f)(z)) = X*(Z(f))- Erdos-Turan's theorem (Theorem I Lip 
then implies that 



A ang (x?(Z(/))<c, 



l lQg / l|£«(/)(*)||sup 



\7~1 & \ \(v,a)\ I ' 

\ D Vnj Res^,...,^(/r, . . . , m— J 

l|£a(/)(s)l|»ui 

n„|R*uy,...,,iS (A™.- -,/»)! ^ 



with c == 2.5619... Lemma E3 implies that log ( - l|Ea(/)(2)l|sup R ^ yr ) is 



bounded above by the quantity 

Il a ll2 lo s k^^t » 

VnjRes^,...,^(/f,...,/-)i^y 

for w = -uAj- G S* n_1 . From the definitions of 8(Z(f)) and of the Erdos-Turan 

|ja|J2 v yJ ' ' 

size rj(f) given in (|2.2p and (|1.3p . respectively, we get 



e(Z(/))<min{l,cvW)}- 

2 

Applying Theorem 12.31 and the fact that the function £5(9 — log(i)) is monotonic 
increasing in the interval (0, 1], we deduce that 

A ang (^(/)) < Mn(j) {9-log(mm{l,cy/^}))^ n - l) mm{l,c^Mf)}-3 

<22n^V§(™- 1 )(18 + log + (7 ? (/)^ 1 ))K"- 1 ) C §r ? (/)^ 

<66n2"(18 + log + (r / (/)- 1 ))l^ 1 )r ? (/)3, 

which gives the bound for the angle discrepancy. For the radius discrepancy, we use 
the bounds given in Theorem ll.l| (|2,2p . and Theorem 12.31 to get, for < e < 1, 

A rad (/,e) < ^2A Tad (xtHZ(f)),e) < — n(f). 

3=1 £ 

This concludes with the proof of the theorem. □ 

We next study a number of basic properties of the Erdos-Turan size. The following 
proposition shows that this notion generalizes the measure of polynomials that appears 
in the statement of Theorem 11.11 



Proposition 3.12. Let d > 1 and f = ao + • • • + a-d xd G C[x] with a^a^ ^ 0. Then 



»,(/) = £ log (J£» 



Proof. The facet resultants of / are 



Res„cn(f° for " = 1 ' 



for v = —1. 



Moreover, D = MV R ([0,d]) = d, D W)1 = 1 for w G S° = {±1}, and |(v,«j)| ■■ = 1 for 
all v, w G {±1}- The formula for r/(/) then boils down to -s log I - l|sup 

\^\a a d \ 

We denote by A n = conv(0, Ci, . . . , e n ) the standard simplex of R n . 
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Proposition 3.13. Let Ai, ■ ■ ■ ,A n be family of non-empty finite subsets of 7U 1 such 
that MV(Qi, . . . , Q n ) > 1 with Qi = conv(„4j). Let fi, ■ ■ ■ , f n G < C[x 1 , . . . , x^ 1 ] with 
supp(/j) C Ai and such that Res^ i ... i ^(/^, . . . , /") / for all t)£Z™\ {0}. 

(1) 7/(/)<00. 

(2) Let Ai,...,A n G C x . Then r?(Ai/i, . . . , A„/„) = r/(/i, . . . ,/ n ). 

(3) £e£ d,- G Z>i and bj G Z n suc/i i/iai N(/,-) C dj-A n + fy, j = 1, . . . , n. Then 



+ E -^ log+ ' Res ^->^> ■ • • > £) _1 i) • 

the second sum being taken over all primitive vectors v G Z n . 1/ /i, . . . , / n G 

Z[x 1 , . . . , x^ 1 ], i/ien 



m< (ti + y^T) ( n"=i rfj) ^ log n/j 

^ JJ - MV(Qi,...,Q„) ^ d,- 



Proof. The statement of ([I]) is clear, since n(f) is defined as the supremum of a 
continuous function over the compact set S" 1-1 . 

For (J2]), let a G Z n \ {0}. As explained in the proof of Theorem ll.4| the product 
of the leading and the constant coefficients of E a (f) is equal to 

±l[Res A1 ,..., A «(f?,...JZf^. 

V 

Hence, the denominator in the definition of n(f) is multihomogeneous in the coeffi- 
cients of each /», of partial degrees equal to 1 1 ct 1 1 2 times those of E a (f). Hence, 

ij— fr de gf 1 (- E «(/)) = MV a-L Ua(Ql), ■ ■ ■ ,^ W {Qj-l),^a(Qj+l), ■ ■ ■ ,7T (Qn)) = A»J 

\\ a \\2 

for w = TTTi-, which implies the statement. 

I! a l|2 

©: let w G S n_1 . Then vr l(; (N(/ i )) C n w (djA n + bj). Due to the monotonicity 
of the mixed volume with respect to the inclusion, plus its properties of homogeneity 
and invariance under translation, we deduce that, for j = 1, . . . , n, 

(3.14) MV w x({^(N(/,))}^) < MV w± ({ir w (d e A n + b t )} w ) 

<(n-l)!(n^)voV(7r w (A™)). 

The projected simplex ^(A n ) can be covered by the union of the projection of its 
facets. One the facets of A n has (n — l)-dimensional volume equal to / _% , while 
the other n facets have (n — l)-dimensional volume equal to (n _ 1 ^ ! . Since the volume 
cannot increase under orthogonal projections, we have that 

(3.15) voU (^( A «))<-^±^. 

In addition, |(t?,ie}| < ||f||2 since w G S n . Then, the first part of the statement 
follows from (|3.14p . (|3.15p and the definition of r/(f). The second part follows from 
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the fact that, if coefficients of the /j's are integers, then the relevant facet resultants 
are nonzero integers and so their modulus are greater or equal to 1. □ 

Let us consider the statement of Proposition I3.13l| 3|) , in the classical dense case 
Qj = djA n for all j. In this situation, the only primitive vectors v to consider are 
v = a, i = l,...,re, and v = eo := — J27=i e i- Given dj > 1, j = 1, . . . , n, we 
write Resd lv ..,d n for the resultant of n homogeneous polynomials in n variables of 
respective degrees d\, . . . ,d n , as defined in [CLO05J . Given a system of polynomials 
fi, ■ ■ ■ i fn £ C[#i, . . . , x n ] with deg(fj) < dj and i = 0, . . . , n, the facet polynomials 
/■f*, . . . , f%* form a system of n polynomials of degrees di,...,d n . In particular, we 
can evaluate Res<2 lv .. l( f n at these polynomials. If we assume that these resultants are 
nonzero, we obtain 

r?(/) < (n + v^)f: l0gl ^ llsUP + 2r i d . (^log + \Re Sdl ,..., d MT,---Jn°r 1 \ 

n v 

+ ^log + |Re Sdl ,... A (/fV..X i r 1 | ■ 

i=l ' 

In particular, if /i, . . . , f n G Z[xi, . . . , x n ], then 

V(f) <(n + Vn) 



Y^ lo g||/jl|sup 

^ d j 



4. Asymptotic equidistribution 

We will apply here the results in the previous sections to study the asymptotic 
distribution of the roots of a sequence of systems of Laurent polynomials over Z and 
of random systems of Laurent polynomials over C 

First, we will consider polynomials over Z. Let Qi, . . . ,Q n C M. n be a family of 
lattice polytopes such that MV(Qi, • • • , Q n ) > 1- For each d > 1 and i = 1, . . . , n, 
consider the finite subset of Z n given by 

(4.1) Ai, d = dQinZ n . 

Proposition 4.2. For d> 1 let f d = (/i <£, . . . , f n ,d) be a family of Laurent polynomi- 
als in Z[xf x , ..., x^ 1 ] such that supp(/ i)d ) C dQ { and Resjq d ,...,AZ d (fi >d , ■■■, fn, d ) ^ ° 

for all d£Z"\ {0}. Then there is a constant c\ > which does not depend on d such 
that 



Er=i lo g|l/i,rfllsu P V (-. . , + / d 



Aang^/J) < d ( E "= ll0g J l/M " SUP ) ' (l + log 



1 . . . . 2 



(n-1) 



,Er=i io gii/i,rfiisu P 

and, /or any < e < 1, 

A C7(f \ *\ <T „ ^i=l & \\ji,d\\sup 
^rad(^(JdJ,ej < Cl —^ • 

Proof. This follows easily from Theorem 11.41 and Proposition |3"TT3"1 J3"|) . □ 

Proof of Theorem \1.7\ Let notations be as in the statement of Theorem ll.71 In partic- 
ular, vz(f d ) 1S the discrete measure associated to Z(f d ), and v is the measure on (C x ) n 
induced by the Haar probability measure on (S 1 ) 71 . 
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We have to show that, for every continuous function with compact support h, 

lim / hduztf ) = / hdv. 

d-+°°J(C><) n d J(C x ) n 

It is enough to prove the statement for the characteristic function hu of the open sets 
of the form 

(4.3) U = {(zt, ...,Zn)€ (C x ) n | nj < \zj\ < r 2 ,j, a, < argfo) < fa for all j}, 

with < r\j < T2,j < oo, rjj ^ 1 and — ir < aj < f3j < n, since any continuous func- 
tion with compact support can be uniformly approximated by a linear combinations 
of these. 

Consider first the case when U n (S' 1 ) n = 0. Due to the conditions imposed on the 
numbers r^j, there exists e > such that U is disjoint with the set {£ £ C n | 1 — e < 
\£j\ < (1 - ey 1 for all j}. Hence, 

f h ha deg(Z(f d )\u) 

/ (CX) „ hud6z ^ = d» MV(Q) - Arad (^' £ )' 

where we set Q = (Qi,...,Qn) and Z{f d )\ v = ^2 ie \ Z (f d )\nu m ^- Proposition 0T2] 
implies that this integral goes to for d — > oo, which proves the statement in this 
case, since J c „ hu dv = 0. 

Consider now the case when U n (S 1 )™ ^ 0. Set U = {z \ a,j < arg(zj) < 
/3j for all j}. Then 

/ hud6 Z (f d )- hudv= h ud5 z{f) -\\^- — l 



J(c x )' 1 V ~~[ 27r / J(c x )" 



We have 

deg(Z(/ d ) ttij g) -A- /3j - o^ 



:)r 



h u ~ 11^— dd z(f d )< ,„ MVfol "11 



2tt 



d" MV(Q) fjj 2vr 



< A ang (/ d ). 



Again, Proposition 14.21 implies that this integral goes to for d — > oo. On the other 
hand, U \ U is a union of a finite number of subsets U\ of the form (|4.3p such that 
U\ n (S 11 ) 71 = for all I. By the previous considerations, f( CX \ n hui&$z{f d ) ~^d and 



so / (CX)n hjj^j d5 z{fd) ^ d 0. Hence 



■a. 



lim / ^ (W z(/d) = J] ^T^ = / ^ di/ = 0, 

which concludes the proof. □ 

We will now consider random systems of Laurent polynomials with complex coeffi- 
cients. To explain and prove our results, we have to consider metrics and measures on 
projective spaces over C. Let N > 1 and consider the standard Riemannian structure 
on C induced by the Euclidean norm || • [|2. Let S 2 = {z £ C | \\z\\2 = 1} be 
the unit sphere with the induced Riemannian structure. The map S 2N_1 —> ¥(C N ) 
given by (jzqi • • • j zn-i) •->• (zq : ■ ■ ■ : zjv-i) gives a principal bundle with fiber S 1 . The 
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Pubini-Study metric on P(C Ar ) is denned as the unique Riemannian structure such 
that this map is a Riemannian submersion, see |KN69j for details. 

The geodesies of P(C Ar ) coincide with the lines. Hence, we can define a distance 
between two points z\ and Z2 as the length of the line segment joining them, and we 
will denote it by distFS^i,^)- However, it will be more convenient to consider the 
distance function dist := sin(distFs)- This function can be computed with the formula 



(4.4) dist( Zl ,z 2 ) l(Zl ' Z2) 



|£l||2|F2||2' 

for any choice of representatives Z{ G C \ {0}, i = 1,2. 

Lemma 4.5. Let Z\,Z2 G P(C ) with Zi G C , i = 1,2, representatives of these 
points, and such that \\z2W2 = 1- Then dist (z\, z 2) < \\z2 — 2i || 2 - 

Proof. We have that 

\\z2-z1W2 = (z2-zi,z 2 -zi) = l+||ii||l-2Re((zi,2 2 )) > 1+ P1II2 - 2|(2i,2 2 )|- 

On the other hand, the formula (j4.4|) gives dist(zi, z 2 ) 2 = 1 — ( n^y ) • Hence, 

||~ ~ i|2 A -.(^ v \2_(u~n \{zi,Z 2 )\ \ 2 

\\Z2- zi\\ 2 - aist{zi,Z2) = [\\Z1W2 [t^— n — > U, 

v lFl||2 J 

which proves the statement. □ 

We will need the following Lojasiewicz inequality for a hypersurface of a complex 
projective space. For a homogeneous polynomial / G C[x] of degree d, and a point 
z G P(C N ), the value 

!/(*)! 

II *\\d 
\\Z\\2 

is well-defined. For a subset E C P(C ), we write dist (2, E) for the distance between z 
and E. 

Lemma 4.6. Let f G C[xo, • • • , aJjv-i] ^ e a homogeneous polynomial of degree d > 

and z £ PfC^). TTien 

|/(z)| ^ / _ 1/(8,)!' 



iM sup ,w) dist( ^ (/))t 

I2 V^ePC" ll^lb / 



Il 2 ll2 \x&(C N ) Il x 'll2 

Proof. Let z, a: G P(C Ar ) such that a; ^ V(f). Let z, £ be representatives of these 
points in the sphere S 1 and set fx(t) = f(z + i£) G C[i]. This is a univariate 
polynomial of degree d with leading coefficient fix). Then, there exist & G C, j = 
1, . . . , d, such that /^ = /(£) \\At — £j) and so 

1/(5)1 =i/ S (o)i= i/(2)i n^-i- 

3 
For each j, we have that z + ^£ G V(f). Using Lemma 14.51 we deduce that 

101 = \\(z + ^jx)-z\\2 >dist(z, z + tjx) >dist(z,V(f)). 
We deduce that 

^ = |/(2)| > \f(x)\dist(z,V(f)) d = \l^dist(z,V(f)) d . 

I y IIU ™ (J/ 

Since this holds for all x ^ V(f), the result follows. □ 
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Let /xfs denote the measure on P(C ) induced by the Fubini-Study metric. Then 



r N-l 



/iFs(P(C )) = (%_-i\\ ■ We will consider the normalized measure given by 



pv=i5 



(N-iy. 



A result of Beltran and Pardo [BP071 Theorem 1] shows that, for a hypersurface 
H £ P(C ) of degree d, the normalized measure of the tube around H of radius p > 
is bounded above by 

15d(N - 1)V- 
Applying this result, we deduce the following bound for the volume of the set where a 
polynomial can take small values. For 5 > and a homogeneous polynomial / £ C[x], 
we consider the subset of P(C Ar ) given by 

(4.7) V(f) s = \zeP(C N ) \IW <S \ 

*• \\ z \\2 

Proposition 4.8. Let 5 > and f £ C[xo, ■ ■ ■ j^jv-i] a homogeneous polynomial of 
degree d > 1. T/ien 

r3 f 5 \ 3 



mCW)*) < 15dA^ d 

MIJ llsup' 

In particular, if f £ Z[xo, • • • , scjv-i]> ^en //(V(/)<0 < 15diV 3 £d. 
Proof. Let 2 £ V(/)<s. Using Lemma 1431 we deduce that 



5 > ( sup l -p&P\ dist(z, V(/)) d > ™^ dist(z, V(f)Y 
\ * \\x\\2 J N2 



Hence, 

dist(z,T/(/)) < m 

' 'I sup, 

and so V(f)s is contained in the tube around V(f) of radius iVa ( TT/ni — ) d ' ^ e first 
part of the result follows then from the Beltran-Pardo bound for the volume of this 
tube. The second part follows from the fact that ||/|| S up > I /I > 1) because of the 
inequality (J3.6J) and the fact that the coefficients of / are integer numbers. □ 



Let us keep notations as above and set A-d = (A\ t d, ■ ■ ■ ,A n ,d) with Ai t d as in (J4.1D , 
Each point of the projective space P(C^ d ) can be identified with a system f d = 
(fi,d) ••■■> fn,d) °f Laurent polynomials such that supp(/j^) £ dQi, i = 1, . . . , n, mod- 
ulo a multiplicative scalar. The associated cycle Z(f d ) is well-defined, since it does 
not depend on this multiplicative scalar. 

Set fid for the normalized Fubini-Study measure on ¥(C^ d ) and let gd- P(C^ d ) — > 
M>o be a probability density function, that is, a //^-measurable function with 

/ g d d/x d = l. 

Let f d be a random system of Laurent polynomials with supped) £ dQi, i = 1, . . . ,n, 
distributed according to the probability law given by g d with respect to pd- We 
can then consider the angle discrepancy of Z(f d ) and, for < e < 1, the radius 
discrepancy of Z(f d ) with respect to e, as random variables on P(C^ d ). We denote 
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by E{A ang {Z{f d ));g d ) and E(A vad (Z(f d ),e);g d ) the expected value of these random 
variables. 

Theorem 4.9. For each d > 1, let g d be a probability density function on P(C d ) and 
fd = (/i,d> ■ ■ ■ > fn,d) a random system of Laurent polynomials with supp(fi >d ) C dQi, 
i = 1, . . . , n, distributed according to the probability law given by gd with respect to /id- 
Assume that the sequence (g d ) d >i is uniformly bounded. Then there is a constant 
C2 > which does not depend on d such that 

los(d+ l)h-h 

E(A ang (Z(f d ));g d ) < c 2 gl > 

d-i 

and, for any < e < 1, 

log(d + l) 



E(A Tad (Z(f d ),e);g d )<c 2 



ed 
In particular, lim d _ >0O E(A ang (Z(/ d ));p d ) = and ]im d ^. 00 E(A iad (Z(f d ),£);g d ) = 0. 

Proof We first estimate the expected value of the angular discrepancy, which is given 
by the formula 

E(A ang (Z(f d ));g d )= [ A ang (Z(f d ))g d (f)d^ d . 

•/P(C*rf) 

Consider the Minkowski sum Q = Yj7=i Qii w hich is a lattice polytope on M. n 
of dimension n because of the assumption that the mixed volume of Qi, . . . ,Q n is 
positive. For each primitive vector d£Z" which is an inner normal to a facet of Q, 
consider the facet resultant 

Rv,d = Res ^r id ,-,^ lid G %[ u i, ... ,u n ], 

where ui is a group of j^Ai^d variables. Proposition I3,13lf 2j) implies that its total 
degree is bounded by deg(R v> d) = c v d n ~ x for a constant c v independent of d. Its total 
number of variables is #«4<2 = X^-^M = Si dQi n Z n . This number can be bounded 
by c%d n for a constant C3 independent of d. 

Set 6 Vtd = d- 2ndc ^ R -' d \ Consider the subset V{R v4 )& v . d C P(C Ad ) as defined 
in (|4.7p and put 

U d = [jV(R Vjd )s Vid , 

V 

the union being over all primitive inner normal vectors to facets of Q. Using the fact 
that < A ang (Z(f d )) < 1 and the hypothesis that the functions g d are uniformly 
bounded, we deduce that 

< / A ang (Z(/ d ))g rf (/)d/i rf < ( sup g d (f d ))fj, d (U d ) < c^ d {U d ) 
Ju d v f d ' 

for a constant C4 independent of d. By Proposition 14.81 

(4.10) fi d (U d ) < Z Z^(y(R V)d ) Svtd ) < ^15deg(i? M )(#Ai) 3 C 

V V 

< 15(5^cd n - 1 )(c 3 d n ) 3 d- 4n = cad.- 1 , 

V 

with c 5 = 15( Y^ v Ct,)c|. Hence, fJ, d (U d ) -> d and so J Ud A axig (Z(f d ))g d (f d ) dfx d -> d 
as well. 



,3 r dc s(R v ,d) 
',d 
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Let f d G P(C' Ad ) \ U d and choose a representative f d = (f lyd , ..., j n4 ) G C Ad \ {0} 
with H/dlb = 1- By Proposition l3.13H[2j) . r](f d ) = rj(f d ). Note that the Minkowski 
sum Yli dQi coincides with dQ. Hence, the only non trivial facet resultants of the 
family of finite sets A\ td , . . . ,A ntd are those of the form R vd as considered above. As 
before, denote with A n the standard simplex of M. n . Choose e > 1 and 6j G Z n such 
that Qi C eA™ + bj for all i. Hence, dA%d C deQi + dbi for all i. Proposition !3.13T fHj) 
implies that there is a constant cq independent of d such that 

1 n 1 

V{fd) ~ d"MV(Q) (( de )" _1 ^ + ^) E lo S HAdlU + 2 E in* lo g + I^C/S)' 

1 n 

< l^TTTTT^ ({de) n -\n + y/n) V log(#A,d) +«V 1Mb deg(i2„, d ) log(d) 



< 



d™ MV(Q) 
log(d + 1) 



= 1 



C6- 



d ' I 

the second and fourth sums being over the primitive inner normals v to the facets 
of Q. Here, we used the fact that ||/t,d||sup < #Ai t d for f d in the unit sphere of 
C"^, the definition of the set U d , the bound #Ai t d < #Ad < c^d n and the inequality 
deg(R v> d) < c„d n_1 that we explained before. 

1 n — 1 

Using Theorem 11.41 and the fact that the function ts log ( j) 3 is increasing for 
small values of t > 0, we deduce that, for f d G P(C ,Ad ) \ U d , 



Hn-l) 



A ang (Z(/,)) < 07i7(/ d )*log (^y)' 

^ /log(d+l)\f. / d \lin-l) log(d+l)§"-3 
- C 9 j lo S i / . , 7\ < Cio 1 

V d J \log(d+l)J dl 

for suitable constants cj, c%, eg and c\$. This proves the first part of the statement. 
For the radius discrepancy, we proceed in a similar way: given e > 0, we write 
E(A ra( j(/ rf ,e); ^) as an integral, which we split into two parts. We bound the first 
using that < A ra( j(/ d ,e) < 1 and the estimate f)4.10[) . while the second integral can 
be bounded by applying Theorem 11.41 □ 

Proof of Theorem \1.8l The proof is similar to the one given for Theorem 11.71 Write 
v d = rf n MVtQ) f° r short, where K(Z(f d );g d ) is the expected zero density measure of 
fd ■ To prove the statement, it is enough to show that, for all subsets U as in (|4.3p . 

d->oo I otherwise. 

If U n (S 1 ) n = 0, then there exists e > such that 

deg(Z(f d )\u) < deg(Z(/ d ))A rad (/„ £ ) < d n MY(Q)A vad (f d ,e) 

and so fd{U) < E(A ra d(/^, e);gd)- Theorem 14.91 then implies that ]hn.d->oo u d(U) = 
= i/(l7). 

If U n (S 1 ) 71 7^ 0, we set U = {z \ otj < arg(zj) < j3j for all j}, and then we have 

"M - n ^ = (»*<& - n h ^ L ) - «& \ ^ 
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Set Rd = Y\ v R Vt d for the product of all non-trivial facet resultants of the family of 
finite sets A\, ■ ■ ■ , A n . Then 



Vd\ 






PCCAiJXV^) 



d™ MV(Q) -"-J- 2vr 



j=l "ir^ «.;\v (ii. di \ »/ •_■£ 



< / a A ang (/ d )5 d (/ d )d// d . 

■McM 



/P^d) 

We have that (7 \ U is a union of a finite number of subsets U\ of the form (|4.3|) such 
that Ui n (5 )™ = for all I. By the previous considerations, lim^oo v d {U{) = and 
so lim^oo v d {U \ U) = 0. Theorem 14.91 then implies that 

n R — 
lim v d (U) = lim ^(17) = TT ^— ^ = *,([/). 

□ 

Remark 4.11. It is not clear to us whether the upper bound in Proposition 14. 81 for the 
volume of the set V(f)s is sharp or not. It would be interesting to clarify this point, 
as a qualititive improvement on this bound might enlarge the range of applicability 
of theorems 14.91 and FO 



Remark 4.12. In some situations, it might be interesting to consider probability 
distributions on the complex linear space C^ d rather than on P(C"^). For a point f d G 
P(C"^), the associated cycle Z(f d ) does not depend on the choice of a representative in 
C d for this point and, a fortiori, the same holds for the angle and radius discrepancies 
of Z(f d ). Hence, one might consider random variables on C d arising from this 
cycle as random variables on P(C d ), by applying Federer's coarea formula (see for 
instance |BP07j Theorem 20]). 

In precise terms, the normal Jacobian of the map w. C^ d \ {0} — > P(C^ d ) with 
respect to the Euclidean structure on C^ d \ {0} and the Fubini-Study one on P(O^) 
is given, for g G C"^ \ {0}, by 

Nj g CT = ii<,ir 27Vd . 

with Nd = #«4d — 1. Given a probability density function Gd'- C^ — > R, one might 
derive a corresponding probability density function on P(C^ d ) by integrating along 
the fibers of w as 



7r Nd r „„ 

(4-13) g d (f d ) = j^ / G d (g)\\g\\l Nd dw-\fd), 

where dw~ 1 (f ' d ) is the volume form of the fiber zu~ 1 (f d )- The probability distribution 
given by Gd of, for instance, the angle discrepancy, can then be computed, for any 
Borel subset I C [0, 1], as 

Prob(A ang (Z(/ d )) G I; G d ) = I g d (f d ) ^ d , 



with Aj^(J) = {f d G F(C Ad ) | Aang(Z(/ d )) G I}. This is a consequence of the 
coarea formula. 
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Example 4.14. Let f d = (/id, . . . , f n d) be a random system of Laurent polynomi- 
als with supp(fi : d) C dQi whose coefficients {fi,d,a\a&A id are independent complex 
Gaussian random variables with mean and variance 1. This is a probability distri- 
bution on C^ d whose density function is defined, for f d € C^ d , as 



G d (fo 



n 



" \Ji,d,a I 



1 



=la6A, 



7T 



#A d 



-Il/ d ||| 



The random cycle Z(f d ) might be described by a probability distribution on P(C^). 
The corresponding density function is the constant function g d = 1. This can be seen 
by computing the integral along the fibers (|4.13p , or simply by observing that Gd is a 
function of the radius H/Jh- 

Theorem 11.81 implies then that the sequence of roots of f d converge weakly to the 
equidistribution on (S 1 ) 71 when d — > 00. In this way, we recover a result of Bloom and 
Shiffman [BS071 Example 3.5]. 
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